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ABSTRACT 

We present quantitative predictions for the detectability of individual Galactic dark matter subhalos 
in gamma-rays from dark matter pair annihilations in their centers. Our method is based on a 
hybrid approach, employing the highest resolution numerical simulations available (including the 
recently completed one billion particle Via Lactca II simulation) as well as analytical models for 
the extrapolation beyond the simulations' resolution limit. We include a self-consistent treatment of 
subhalo boost factors, motivated by our numerical results, and a realistic treatment of the expected 
backgrounds that individual subhalos must outshine. We show that for reasonable values of the 



dark matter particle physics parameters (M^ ~ 50 — 500 GeV and 



10- 



GLAST may very well discover a few, even up to several dozen, such subhalos, at 5 cr significance, 
and some at more than 20 cr. We predict that the majority of luminous sources would be resolved 
with GLAST's expected angular resolution. For most observer locations the angular distribution of 
detectable subhalos is consistent with a uniform distribution across the sky. The brightest subhalos 
tend to be massive (median V^nax of 24 kms^^) and therefore likely hosts of dwarf galaxies, but many 
subhalos with Umax as low as 5 kms~^ are also visible. Typically detectable subhalos are 20 — 40 kpc 
from the observer, and only a small fraction are closer than 10 kpc. The total number of observable 
subhalos has not yet converged in our simulations, and we estimate that we may be missing up to 3/4 
of all detectable subhalos. 

Subject headings: dark matter - Galaxy: structure - galaxies: halos - gamma rays: theory - methods: 
n-body simulations 



> 

o 
oo 
O 



- 

X 



1. INTRODUCTION 

Revealing the nature of dark matter is fundamental 
to cosmology and particle physics. In the standard cos- 
mological paradigm of structure formation ( ACDM) , the 
universe is dominated by cold, coUisionless dark matter 
(CDM), and endowed with initial density perturbations 
via quantum fluctuations during inflation. In this model 
galaxies form hierarchically, with low-mass objects ("ha- 
los") collapsing earlier and merging to form larger and 
larger systems over time. Small halos collapse at high 
redshift when the universe is very dense, so their central 
densities are correspondingly high. When these halos 
merge into larger hosts, their high densities allow them 
to resist the strong tidal forces that act to destroy them. 
It is therefore a clear, unique prediction of ACDM that 
galaxies are embedded in massive, extended dark matter 
halos teeming with self-bound substructure or "subha- 
los". 

If the dark matter is in the form of a relic particle 
once in thermal equilibrium in the early universe, like the 
neutralino in SUSY, then substructures will be lit up by 
pair annihilations of these particles into standard model 
particles, whose subsequent decay results in gamma-ray 
emissions. The signal depends on unknown conjectured 
particle physics properties and poorly known astrophysi- 
cal parameters. Particle physics uncertainties include the 
type of particle (axion, neutralino, Kaluza-Klein particle, 
etc.), its mass, and its pair annihilation cross-section. 
From the astrophysical side, since the annihilation rate 
is proportional to density squared, the predicted flux de- 
pends sensitively on the dumpiness of the mass distribu- 



tion. 

The detection of annihilation radiation from DM 
clumps is surely one of the most exciting discoveries 
that the upcoming Gamma-ray Large Area Space Tele- 
scope (GLAST) could make. Currently scheduled to be 
launched in May of 2008, GLAST is NASA's successor to 
the Compton Gamma Ray Observatory, whose EGRET 
instrument conducted the first all-sky gamma-ray sur- 
vey and disc overed 271 sources, 170 of which remain 
unidentified (jHartman et al.lll999f ). GLAST will carry 
two instruments, the GLAST Burst Monitor (GBM), de- 
signed to detect flashes from gamma-ray bursts and solar 
flares, and the Large Area Telescope (LAT), the main 
survey instrument. The LAT consists of several layers of 
high-precision silicon tracking detectors for determining 
the direction of an incident gamma-ray, a cesium-iodide 
calorimeter to measure its total energy, and an anticoin- 
cidence detector for cosmic ray rejection. 

With a field of view of about 2.4 sr and a 90 minute 
orbital period, the LAT will survey the entire sky daily. 
Its peak effective area exceeds 8000 cm^ beyond IGeV 
and its angular resolution ranges from < 3.5° at 100 
MeV down to < 9 arcmin above 10 GeV. Compared to 
EGRET, the LAT has a 4-5 times larger field of view, 
more than 5 times larger peak effective area, 5-40 times 
higher angular resolution above 1 GeV, and ~ 30 times 
better point source sensitivity. Additionally, the LAT 
detector is sensitive to gamma-rays out to 300 GeV, ten 
times higher than EGRET's limit. Given these improve- 
ments it is not unreasonable to hope that GLAST will 
discover previously unknown sources, one of which may 
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TABLE 1 
Simulation Parameters 



Name 


^box 


e 


Zi 


-^hircs 




»'200 


M200 


Vniax 


''"Vmax 


A^sub 




(Mpc) 


(pc) 






(Mo) 


(kpc) 


(Mq) 


(km s^-"-) 


(kpc) 




VL-I 


90.0 


90.0 


48.4 


2.34 X 10** 


2.1 X 10"* 


389 


1.77 X 10^2 


181 


69 


9,224 


VL-II 


40.0 


40.0 


104.3 


1.09 X 10^ 


4.1 X 10^ 


402 


1.93 X 10^2 


201 


60 


53,653 



Note. — Box size Lboxi (spline) softening length e, initial redshift Zi, number N-^^^^^g and mass Afhircs of high resolution dark matter 
particles, host halo r200, and Vmax, and number of subhalos within r20O A^sub for the VL-I and VL-II simulations. Force softening 

lengths e are constant in physical units back to 2 = 9 and constant in comoving units before. 



be annihilating DM in the centers of Galactic subhalos. 

What region in the sky is mostly likely to allow for 
a detection of gamma-rays from DM annihilations is a 
hotly contested question. The Galactic center (GC) is 
likely the largest nearby DM overdensity and a number 
of studies have advocated it as the best place to look for 
a signal (^erezin skv ct al. 1994; Berg strom ct al. 199^ 
[200I lUllio eTan 120021: ICesarini et allllool . A major 
drawback for this scenario, however, is the large astro- 
physical gamma-ray background near the GC. Further- 
more, the lack of a spectral break below ~ 30 TeV in 
the recent H.E.S.S (|Aharoman et al.|[2006l ) and MAGIC 
(| Albert et al.ll2006[ ) measurements of the very high en- 
ergy gamma-ray emission from the GC severely limits 
the contribution that annihilating DM particles with 
masses l ess than 0(10 TeV ) can make to this signal. 
Instead, iStoehr et al.l ()2003f ) suggested to focus obser- 
vations on an annulus between 25° and 35° from the 
Galactic center, since this would maximize the chances 
of detection by reducing the background while simul- 
taneously avoiding numerical uncertainties in the cen- 
tral DM density profile. Perhaps more promising is the 
possibility of detecting DM annihilation from the cen- 
ters of Galactic subhalos, either from one of the known 
Milky Way dwarf satellite galaxies or from a nearby 
dark clump (iBer gstrom et al.l Il999l: iBaltz et all l2000l: 
Calcaneo-R oldan fc Moorg 20001 iTasitsiomi fc Olinto 



2002| :IStoehr et al.ll2 003VTavk3r fc Silk"2003': 'Eva ns et al 



20041: lAloi sio et al. 2 004: Koushiappas et all 12004 



Koushiapp as 2006; Die mand et al.H2007at IS~ gari et al.l 



20071j : iPieri et al.i i2008h . Lastly, an additional source 



population might arise from local DM density enhance- 
ments (mini-spikes) aro und black hole remna nts of the 
first generation of stars (jBertone et al.l[2"005bl ). 

If GLAST detects photons originating in pair annihi- 
lations of DM particles in the centers of subhalos, will it 
be possible to distinguish this signal from conventional 
astrophysical gamma-ray sources? While a definite an- 
swer to this quest i on wi ll have to await actual GLAST 
data, iBaltz et al] (j2007[ ) have argued that this should 
be possible based on four criteria: (i) a hadronic spec- 
trum from monochromatic quarks, (ii) lack of time vari- 
ability, (iii) spatial extent, and (iv) a lack of emission 
at shorter wavelengths, except for very diffuse inverse 
Compton and synchrotron radi ation (jBaltz fc Waill2003 : 
IColafrancesco et al.|[200ll2007[) . Gamma -ray pulsars are 
probably the most problematic astrophysical sources in 
this regard, owing to the similarity of their spectra to 
a typical DM annihilation spectrum. Fortunately they 
tend to lie in the Galactic plane, often exhibit X-ray 
counterparts, and are point-like. 

The observability of DM annihilation radiation orig- 
inating in Galactic subhalos depends on their abun- 



dance, distribution, and internal structure, properties 
which must be determined by numerical simulation, since 
they are the result of structure formation in the highly 
non-linear regime. Previous investigations have typi- 
cally only indirectly made use of numerical simulations, 
employing (semi-)analytic models that have been cali- 
brated to published numerical results. Often these re- 
sults were derived from simulations that do not resolve 
the relevant sub-galactic scales, and it is not clear that 
a simple extrapolation is justified, since not all substruc- 
ture p roperties are scale invariant . Some pas t work 
(ICalckn co-Roldan & Moore 20 00l: IStoehr et all l2003l : 
iDiemand'et all l2006l : lAthanassoula et al.l l2008f l has di- 
rectly used numerical simulations, but these studies ei- 
ther suffered from insufficient resolution, didn't correct 
for effects from the population of unresolved subhalos, 
or didn't realistically account for the gamma-ray back- 
grounds against which individual subhalos must com- 
pete. 

Ground based gamma-ray detectors, i.e. atmospheric 
Ccrenkov telescopes (ACT) like H.E.S.S., VERITAS, 
MAGIC, and STACEE will provide complementary in- 
formation to GLAST's observations. However, owing to 
the comparatively small field of view of ACT's, the detec- 
tion of individual subhalos with such observations would 
either have to rely on serendipity to provide a detectable 
subhalo in the surveyed portion of the sky, or they would 
have to target a known likely source such as a nearby 
dark matter dominated dwarf ga l axv (e.g. lAlbert et al.] 
[2008t iSanchez-Conde et al.l [2007t iDriscoUet al.l l2007f) . 
As advocated by Koushiappas et aP (|2004f ). a better 
strategy might be to locate sources with GLAST and 
then conduct follow-up observations with an ACT to 
measure the gamma-ray spectrum out to higher energies. 

Motivated by the imminent launch of the GLAST 
satellite, we present here a comprehensive analysis of 
the detectability of individual DM subhalos, taylored 
to GLAST expectations. Our method is based on a 
hybrid approach, combining the highest resolution nu- 
merical simulations of the distribution of Galactic DM 
substructure (described in § 12. ip with analytical mod- 
els for the extrapolation of the substructure hierarchy 
below the simulations' resolution limit. We include a 
self-consistent treatment of subhalo boost factors (i) l2.2p . 
motivated by our numerical results, and a realistic treat- 
ment of the expected backgrounds that individual subha- 
los must outshine (? 13. ip . We consider a physically mo- 
tivated region of particle physics parameter space (§ 13. 2p 
and check that our results don't violate existing EGRET 
constraints 13. 3p . Our analysis results in quantitative 
predictions for the number of observable subhalos as a 
function of particle physics parameters (§ 13. 4p . and we 
also consider the possibility of detecting microhalos be- 
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low 1 Mq (§[33]). A discussion of our results and conclu- 
sions can be found in § |4l 

2. MODELS AND METHODOLOGY 

The number of DM annihilation gamma-ray photons 
from a solid angle Afi along a given line of sight {9, 0) 
over an integration time of Tgxp is given by 
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'2M2 



dE 



A,s{E)dE 



X / piiydi. 



(1) 



los 



This expression contains terms that depend a) on the 
properties of the detector (the angular resolution Ail, 
the energy dependent effective area Aefi{E), the energy 
threshold Eth, and the total exposure time Toxp), b) on 
the chosen particle physics model describing the nature 
of the DM particle (i.e. its mass and thermally av- 
eraged velocity- weighted annihilation cross section {av) , 
and the photon spectrum due to a single annihilation 
event dN^/dE), and c) on the spatial distribution of the 
DM within our halo. In this section we will focus on 
the astrophysical contribution, while incorporating par- 
ticle physics and detector properties for the predictions 
of subhalo detectability in the following section. 



presented in this paper have been performed for both 
simulations, but the final results are based on the higher 
resolution VL-II. 

Given the position of a fiducial observer, located 8 kpc 
from the host halo center, we bin up the sky into a grid of 
2400 X 1200 pixels, equally spaced in longitude ((/>) and 
in the cosine of the co-latitude {cos 9). This pixel size 
was chosen to approximate the best angular resolution 
of GLAST's LAT detector, equal to about 9 arcmin (sec 
Section [3]) . Each pixel corresponds to a solid angle of 
An = (27r)/2400 x 2/1200 = 4.363 x 10"^ sr. We then 
calculate each particle's angular coordinates on the sky 
{(j)i,cos6i), and sum up the fluxes, Fi = mpPi/Aird^, of 
all particles in a given pixel. Here di is the distance from 
the i**^ particle to the observer, pi is its density, measured 
over a 32 particle SPH-kernel, and we have approximated 
/ pfdV as rUpPi, which is appropriate for a collection of 
discrete simulated DM particles of mass rrip. In order 
to minimize shot noise arising from particle discreteness 
we smooth each particle's contribution using a projected 
SPH- kernel of angular width hi/di, where hi is half the 
radius encompassing all 32 neighboring particles. The 
resulting map is presented in the top left panel of Fig. [1] 
in units of GeV^ kpc cm^^ sr^^ and depends only the 
choice of the observer's location and the dark matter dis- 
tribution. The flux from only those particles belonging 
to subhalos is shown in the top right panel. 



2.1. Simulations 

The "Via Lactea" simulations follow the formation and 
evolution of the dark matter substructure of a Milky- 
Way-scale halo in a cosmological volume. The first of 
these simulations (VL-I) consists of 234 million high 
resolution particles, covering the virial volume and its 
immediate surroundings of a M200 = 1-77 x 10^^ Mq 
halo. With a particle mass of ~ 21,OOOM0, this sim- 
ulation resolves around 10,000 subhalos within r2oo = 
388 kpc^. The global z = properties of the host 
halo and the substruct ure population were presented in 
iDiemand et all ()2007al ). the joint temporal evolution of 
host halo and s ubstruc ture properties was discussed in 
IDiemand et al.l ()2007H ). and an analysis of the shapes, 
alignment, and spatial distribu tion of the subhalos can 
be found in lKuhlen et al.l (|2007f ). Recently we have com- 
pleted "Via Lactea 11", the next genera tion of this se- 
ries of simulations (jDiemand et al.lT2008D . At a cost of 
over 1 million cpu-hours on Oak Ridge National Lab's 
Jaguar supercomputer, VL-II employs just over one bil- 
lion 4, 100 M0 particles to model the formation of a 
-^^200 = 1-93 X 10^2 Mq halo and its substructure. For 
this simulation we use d for the first tim e a novel adap- 
tive time-step method (jZemp et al.ll2007() . which assigns 
time-steps equal to 1/16 of the local dynamical time (but 
no shorter than 268,000 yr). This allows us to accurately 
resolve the density structure much farther into the cen- 
tral regions of the host halo. With VL-II we are able 
to resolve more than 50,000 individual subhalos within 
the host halo's r2oo = 402 kpc. Roughly 2,000 of these 
arc found within 50 kpc of the center, and 20 within the 
solar circle. For reference the parameters of these two 
simulations are summarized in Table [T] The calculations 

^ Note that wc define r200 as the radius enclosing 200 times the 
mean density of the universe 



2.2. Subhalo Boost Factor 

The hierarchical nature of CDM structure formation 
implies that substructure should be expected not only 
in the host halo, but also in individual subhalos. In- 
deed, in VL-II we are able to resolve some of these 
sub-subhalos in the most massive of our subhalos. In 
IDiemand etlHI (|2008[ ) we show that the mean abundance 
of sub-substructure is consistent with a scaled-down ver- 
sion of the VL-II host halo, i.e. once the sub-subhalos' 
Knax's are scaled by the Vmax of their host subhalo, the 
cumulative Knax-function (iV(> Knax)) agrees with the 
overall subhalo ymax-function. A density squared pro- 
jection of four of the most massive subhalos is shown in 
Fig. O Only particles within the subhalo's outer radius 
Vsuh have been included in the projection. ^ The subha- 
los' riooo: the radius at which the mean enclosed density 
is equal to 1000 times the mean matter density of the 
universe, is overplotted with a dashed line. 

Even with VL-II's impressive dynamical range of ^ 10^ 
in mass, however, we still resolve only a small fraction of 
the total expected DM substructure hierarchy. In prin- 
ciple this hierarchy could extend all the way down to 
the cut-off in the matter power spectrum, set by col- 
lision al damping and free streaming in the early uni- 
verse (|Green et al.ll2005HLoeb fc Zaldarriagall2005f) . For 
WIMP dark matter, typical kinetic decoupling tempera- 
tures range from MeV to GeV, co rresponding to cut-o ff 
masses of toq = 10"^^ to 10"'' Mg (jProfumo et al.ll2006f) . 
some 10 to 20 orders of magnitude below VL-II's mass 
resolution. Since the annihilation rate goes as density 
squared, any dumpiness will lead to an enhancement of 
the total luminosity compared to a smooth mass distri- 
bution. 

^ '"sub is defined here as the radius at which the total density 
becomes roughly constant. 
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Fig. 1. — AUsky maps of the annihilation signal from VL-II, for an observer 8 kpc from the halo center along the host halo's intermediate 
principal axis (top two rows) and along its major axis (bottom row). The maps in the left panel show the total annihilation signal from all 
DM particles within r200i ™ the right panel only the signal from subhalo particles. In the top row we show the uncorrected signal directly 
from the simulation, in the bottom two rows the halo center, indicated with a small white ellipse, has been replaced with an artificial 
p oc cusp, and a mass-dependent boost factor (for mo = 10~^ Mq, a. = 2.0) has been applied to the subhalos. 



The unresolved portion of the substructure hierarchy 
will thus lead to a boost in the true annihilation lumi- 
nosity of individual subhalos, compared to the value de- 
termined directly from the numerical simulation. The 
magnitude of this boost will depend sensitively on the 
properties of the subhalo population, in particular on the 
slope and low mass cut-off of the subhalo mass function 
and on the subhalo concentration-mass relation. Our an- 
alytic boost factor mo del is based on the one presented in 
IStrigari et al.l (|2007af l. The true luminosity of a subhalo 
of mass M is given by 

L(il/) = [1 + B{M)\ L{M), (2) 



where 



is the luminosity of the subhalo determined from all sim- 
ulation particles within the subhalo's outer radius and pi 
is the density of the i^^ such particle. A subhalo's B{M) 
can be calculated by integrating luminosities over its own 
sub-subhalo population: 

B(M)^-. / —L(m)dm (4) 

^ ' L{M) J,no dm ^ ' 
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= 218 kpc 




= 21,4 km s'"' 


^sub 


= 46,8 kpt: 
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Fig. 2. — Sub-substructure in four of VL-II's most massive subhalos. Shown are projections of for all particles within a subhalo's 
outer radius r^^i^. The dashed circle indicates the subhalo's riooo- The clumpy sub-substructure boosts the total annihilation luminosity 
of its host subhalo. 



L{M) 



dN 
dm 



[1 + B{m)]L{m)dm. (5) 



Here dN/dm is the sub-subhalo mass function, and the 
integration extends from mo, the low mass cut-off of 
the substructure hierarchy, to an upper hmit of mi = 
minjlO^ Mq, O.IM}, such that only substructure below 
VL-II's resolution limit of ~ 10^ Mq contribute. For 
subhalos below 10^ Mq we cap the integration at Q.IM 
under the assumption that efficient dynamical friction 
would have lead to the tidal destruction of larger sub- 
subhalos. For a power law substructure mass function 



dN/dm = A/M{m/M)-°', Eq. becomes 



B{M) 



A 



In nil 



/ m \ - 

\m) 



[1 + _B(m,)] i(m)(ilnm. 



L{M) Jin mo 

Motivated b y ou r num erical simulations 
(jPiemand et al.l l2004al l2007af) and semi-analytic 
studies (jZentner fc BuUockl l2003f ). we normalize the 
sub-subhalo mass function by setting the mass fraction 
in subclumps with masses 10~^ < m/M < 10~^ equal 
to 10%. 

For the determination of L{M) we have assumed an 
NFW density profile, in which case the total annihilation 
luminosity of a halo of mass M and concentration c = 
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Fig. 3. — The annihilation luminosity boost factor due to sub- 
structure below VL-IFs resolution limit versus subhalo mass, for 
different subhalo mass functions. Top panel: Dependence on the 
cutoff mass mo for slope a = 2.0. Bottom panel: Dependence on 
a for mo = 10~^ Mq. 



r 

/ rs is given by 



L{M,c) cx piri cx M 



(7) 



where /(c) = ln( l + c) — c/(l + c). We use the 
iBuUock et a"n(|200ll ) concentration-mass relation for field 
halos, albeit with a somewhat smaller val ue of the nor- 
malization, K = 3.75 ( as suggested by iKuhlen et al.1 
l2005t iMaccio et all 120071 1. For the cosmology used in 
the VL simulations and halos masses between 10^ and 
10^° M0, the c(M) relation is approximately c(M) « 
18(M/1O8M0)-O■o^ which corresponds to L(M) oc 
M°-®^, i.e. the annihilation luminosity scales almost 
linearly with mass, i n agreement with r esults from nu- 
merica l simulations (jStoehr et al.l l2003t iDiemand et al.l 
[20073). Note that in our numerical simulations we find 
systematically higher subhalo concentrations closer to 
the host halo center. This trend docs not affect the mag- 



nitude of the boost factor, but translates to a radial trend 
in subhalo luminosity (see Section 3.1). 

Eq. [5] is solved numerically using the boundary con- 
dition B{ma) = 0. The resulting relation is plotted in 
Fig. m for a = 2.0 and different values of mo in the top 
panel, and for toq = 10^^ Mq and different values of a 
in the bottom panel. Overall we find relatively modest 
boost factors on the order of a few, ranging up to ~ 10 
for the most massive subhalos. Generally more massive 
halos have larger boost factors, simply because their sub- 
halo population covers more of the total subhalo hierar- 
chy. For the same reason, smaller values of toq lead to 
larger boost factors. For a < 2.0 B{M) has a weaker 
mass dependence and is less sensitive to mo, since in this 
case more massive halos are relatively more important. 
Our r esults are in ag r eemen t with the analytic upper lim- 
its ofJStriga^^FayjISQOT^) and the recent calculations 
of lLavalle et all (|2008l)~ 

A fit to the cumulative subhalo mass function in our 
simula tions is consistent with a = 2 (jPiemand et al.l 
l2007a() , which implies equal mass in subhalos per decade 
of subhalo mass. However, fits to the differential mass 
function tend to favor slightly shallowe r slopes of 1.8—1.9 
(jStoehr et al .1120031 : [Madau et al.ll2008D . possibly because 
they are more sensitive to the lower mass end, where res- 
olution effects may artificially flatten the slope. In this 
work we use a = 2.0 and mo = 10^ Mq as our fiducial 
model, but present results for a range of different a and 
mo. 

2.3. Central Flux Corrections 

The host halo center is another area where our simu- 
lation must be corrected to account for the artificially 
low density caused by the finite numerical resolution 
(IDiemand et al.ll2004bl) Based on numerical convergence 
studies ( Diemand et al.|[2005al) we believe that we can 
trust the radial density profile of the VL-I host halo 
down to Tconv = 3.4x 10~^r2oo = 1-3 kpc (jPiemand et al.l 



l2007a|) . corresponding to about 10° from the center. The 
higher mass resolution and improved time-step criterion 
in VL-II results in a much smaller convergence radius of 
fconv = 380 pc. The flux derived directly from the simu- 
lated particles in VL-II will thus only underestimate the 
true annihilation flux within the inner ^ 2° from the 
center. An additional uncertainty arises from the fact 
that our purely coUisionless DM simulation completely 
neglect the effect of baryons. While this is not a prob- 
lem for the signal from individual subhalos, which are 
small enough that baryonic effects are likely negligible, 
the central region of our host halo most likely would have 
been affected by gas cooling, star formation, and stellar 
dynamical processes. It is not immediately obvious how 
such baryonic effects would alte r the central DM distri- 
bution. Adiabatic c ontraction (jBlumenthal et anil986l : 
iGnedin et al.l[2004a^ would lead to a steepening of the 
central DM density proflle at scales of a few kpc and be- 
low. A recent study of scaling relations in spiral galaxies, 
however, seems to favor models of spiral galaxy formation 
without adiabatic contraction, and suggests that clumpy 
gas accretion might have reduced central DM densities 
(jDutton et al.ll200"7l ). Stirring by a ste llar bar could also 
eject DM from the central regions (jWeinberg fc Kat"3 
l2007l and references therein). On much smaller scales 
(central few pc), the presence of a supermassive black 
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Fig. 4. — The gamma-ray intensity from DM annihilations in 
the smooth host halo versus the angle from the Galactic center, 
directly from the simulated VL-II particles {thick solid line), and 
with an artificial central cusp with central slope equal to 1.0 {dotted 
line) and 1.2 {dashed line). For comparison the smooth host halo 
flux from the VL-I simulation is overplotted {lower thin solid line). 



hole ( SMBH) would lead to the c reation of an r mini- 
cusp ()Gnedin fc Primack|[2004bD . A SMBH binary, on 
the oth er hand, could hav e removed DM from the very 
center (iMerritt et al.ll2002[ ). 

Here we have attempted to correct for the artificially 
low central density by excising all simulated particles 
within Tconv and replacing them with an artificial cusp 
consisting of about 25 million 5 M© particles, distributed 
in an ellipsoid matched to the shape and density of the 



VL-H host halo at r, 



and extended inwards with a 



power law density profile. The resulting radial flux pro- 
file (excluding subhalos) is plotted in Fig. [H for the un- 
corrected case and for central density slopes of -1 and 
-1.2. 

Note that a similar correction should be applied to 
the centers of all the subhalos as well. This, however, 
is computationally very expensive, and we omit a direct 
correction of the particle distribution in favor of applying 
an a posteriori correction to the central pixel of each 
subhalo in the final allsky map. For each subhalo we 
estimate the central surface brightness according to 



(8) 



where J^haio is the flux from the central region subtended 
by solid angle AfJ of a subhalo with NFW sc ale de nsity 
Ps and scale radius rg at a distance D (see Eq. lAlip . We 
estimate ps and from the measured values of T^uax and 
?'Vmax according to the relations 



rvn 



Ps 



2.163 
4.625 

AttG 



(9) 
(10) 



which hold for an NFW density profile. If the central 
pixel of a subhalo has $ < $c, we set it equal to $c- 
Note that we first apply the boost factor correction to 
all particles of a given subhalo, and then apply this flux 



correction to the central pixel. We only correct the cen- 
tral pixel, so subhalos for which subtends more than 
one pixel (1.7% of all VL-II subhalos) will not be fully 
corrected. 

2.4. Corrected Allsky Maps 

In Fig. [T] we show both unmodified and corrected DM 
annihilation flux allsky maps. While the topmost row 
shows the signal as calculated directly from the numer- 
ical simulation, in the center row the subhalo particle 
fluxes have been boosted by B{M), assuming (a, mo) = 
(2.0, 10~^), the central pixel flux correction has been ap- 
plied, and the central region of the host halo has been 
replaced with an artificial density cusp of slope -1. In 
these two rows the observer is located at 8 kpc along the 
intermediate axis of the host halo ellipsoid. In the pan- 
els on the right hand side, we show only the fluxes from 
particles belonging to subhalos, i.e. within a given sub- 
halo's outer radius. The addition of a substructure boost 
factor lifts the brightness of a number of subhalos over 
the level of the diffuse host halo flux. The bottom panel 
shows the corrected case for an observer placed at 8 kpc 
along the major axis of the host halo. In this case the 
observer is deeper inside the host halo mass distribution, 
and the resulting diffuse host halo flux dominates over 
most of the subhalos. 

3. SUBHALO DETECTABILITY 

In order to make predictions for the number of de- 
tectable subhalos, we must convert the simulation fluxes 
presented in Section [5] into actual gamma-ray photon 
fluxes and compare these to th e expected backgro und 
signal. Following previous work (jStoehr et al.ll2003l ). we 
define for each subhalo in our allsky map a "detection 
significance" 

S=^iM^, (11) 

where Ns^i and N^^i are the number of source and back- 
ground gamma-rays received on a given pixel, and the 

1 /2 

sums are over all pixels with signal-to-noise {Ns.i/Nj^ ^ ) 
greater than unity covered by a subhalo. For this pur- 
pose we must make assumptions about the observation 
(the detector's effective area and angular resolution as a 
function of energy, the exposure time and energy thresh- 
old), the astrophysical backgrounds (extragalactic and 
Galactic), and, most importantly, the nature of the DM 
particle (its mass and annihilation cross section, the an- 
nihilation spectrum). 

Here we adopt parameters appropriate for a GLAST 
observation with the Large Area Telescope (LAT): an 
exposure time of Tcxp = 2 years, the expected energy- 
dependent effective area Acff{E), and an energy thresh- 
old of ii^th = 3.0 GcV, above which the flux- weighted 
mean angular resolution 

/^^ Ae{E) Aeff (S) [dN^/dE] dE 



(12) 



Aeft-(i?) [dN^ldE) dE 

is equal to 9 arcmin. The functional form of Acs{E) and 
M{E) was obtained from the GLAST LAT Performance 
webpage'^. 

^ http:/ /www-glast. slac.stanford.edu/softwarc/IS/giast_lat_ 
performance.htm 
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3.1. Backgrounds 

The annihilation signal from individual extended sub- 
halos must compete with a number of diffuse gamma-ray 
backgrounds, of both astrophysical and particle physics 
origin. The Energetic Gamma Ray Experiment Tele- 
scope (EGRET) aboard the Compton Gamma Ray Ob- 
servatory satellite conducted a gamma-ray allsky sur 
vey in the 1990's, detecting a diffuse background con 



sisting of an isotropic extragalactic (ISreekumar et al 



1998 ) as well as a Galactic component (jHunter et al 



19971) . The extragalactic co mponent is probably dom- 



inated by unresolve d blazars (jStecker fc Salamonlll996l : 
iGiommi et al.|[2006l ) , although star-bursting galaxies and 
galaxy clusters might add a substantial contribution 
(|Chiang fc Mukherjed Il998t iDermed l2007l ) . The spec- 
trum of this extragalactic background is well described 
by a power-law ph oton spectrum in the 3 MeV to 100 
GeV energy range (jSreekumar et al.lll998l ): 



: (7.32 ± 0.34) X 10"^ MeV^cm-^s-^sr-i 



E 



451MeV 



(-2.10±0.03) 



(13) 



The diffuse Galactic gamma-ray background is pro- 
duced by energetic interactions of cosmic ray nucleons 
with interstellar medium (ISM) atoms through neutral 
pion production, as well as from energetic electrons by 
inverse Compton upscattering of lower energy photons 
and bremsstrahlung. Modeling this background under 
the assumption that the locally measured electron and 
proton energy spectra are representative of the Galaxy as 
a whole ( "conventional model" ) underpredicts the flux in 
the measured EGRET spectrum a b ove 1 GeV by about a 
factor of two (|Hunter et al.lll997| ). IStrong et all ^OM) 
provide an explanation for this "GeV excess" by model- 
ing cosmic ray propagation within the Galaxy, allowing 
for mild departures in the electron and nucleon spec- 
tra from the conventional model, and using a cosmic 
ray source distribution based on the observed distribu- 
tion of pulsars and supernova remnants i n conjunction 
with a variable CO-to-H2 conversion factor ([Strong et al.l 
l2004bf ). This "optimized model" emphasizes the impor- 
tance of inverse Compton emission and is able to repro- 
duce the EGRET observations in all directions of the 
sky. Here we have employed the GALPROP^ (v50p) cos- 
mic ray propagation co de (jStrong fc Moskalenkol 119981 : 
iMoskalenko et aP l2002l . and references therein) to cal- 
culate an "optimized model" allsky map of the gamma 
ray emissivity due to cosmic ray interactions with the 
Galactic ISM. We have limited the pixel size of this map 
to 0.5°, corresponding to the angular resolution of the 
HI and CO surveys used in GALPROP. 

In addition to these astrophysical backgrounds, indi- 
vidual detectable subhalos must outshine DM annihila- 
tions from the smooth host halo as well as the back- 
ground fro m the popu l ation of individually undetectable 
subhalos (|Pieri et al.l l2008[ ). Note that these unde- 
tectable subhalos do not simply uniformly boost the flux 
from the smooth host halo. Tidal disruption of satellites 
is more effective close to the host halo center, leading to 
a spatial distribution of subhalos that is antibiased with 

* available from http://galprop.stanford.edu/web_galprop/ 
galprop_home . html 




50 100 
ip [degrees] 



150 



Fig. 5. — Diffuse flux due to undetectable subhalos as a func- 
tion of angle ip from the Galactic center, for a number of differ- 
ent subhalo mass functions. The thick lines show models with 
an anti-biased radial distribution, concentrations increasing to- 
wards the host center, and different values of the mass function 
slope a and low mass cutoff mo: (a, mo/MQ) = (2.0, 10~^) 
(solid), (2.0,10-12) (dotted), (2.0,1) (dashed), (1.9,10"'^) (dot- 
dashed), (1.8, 10 ~^) (triple-dot-das hed). The thin solid line repre- 
sents the orig inal lPieri etal] l200g) model (Brof,zo)> with a = 2.0, 
mo = 10~^Mq, an un-biased radial distribution, and no radial 
concentration dependence. The flux from the smooth host halo is 
overplotted with the grey line, see Fig. l4l 



respect to the host halo mass distribution (jKuhlen et al.l 
l2007tlMadau et al.|[2008[) . The resulting background wiU 
have a shallower angular dependence than the smooth 
host halo signal, and can dominate it at large angular 
distances from the center. 

To determine the magnitude and angular dependence 
of this background w e repeat the calculation presented 
in iPieri et al.l ()2008f ). with three important differences. 
Firstly, we use an antibiased subhalo sp atial distribu- 
tion for which nsub(?')/phost(?') oc r (see iKuhlen et al.l 
|2007| ). as opposed to one that follows the host halo mass 
distribution down to some hard cut-off rniin(Af) as in 
IPieri et all l|2008f) . Secondly, we allow a range of values 
of the subhalo mass function slope a and cutoff mass toq. 
This can make a big difference, since by number the pop- 
ulation of individually undetected subhalos is dominated 
by objects with masses close to ruQ. Lastly, we include 
a radial dependence in the subhalo mass-concentration 
relation, motivated by numerical simulations which tend 
to find higher co ncentrations for subhalos clo ser to the 
host halo center (jDiemand et al.ll2007bl . |2008[ ). 



sub 



(Af,i?) 



ROl 
Co 



\r200 



-0.286 



(14) 



where Cg''^(M) is the medi an concentrati o n of subhalo 
of mass M, as given by the iBuUock et al] ()200lD model 
for field halos. With this scaling, subhalos at Rq are 
three times as concentrated as field halos. We also 
include a log-normal s catter around this in edian, with 
width cTiogjoc = 0.14 (jWechsler et al.ll2002[ ). In Fig. [5] 
we present the resulting background flux as a function of 
angle from the halo center, and also show the effe cts of 
our modiflcations on the original iPieri "eTaLl (l2008h pre- 
scription (using their Srof.zo model). A more detailed 
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TABLE 2 

SuBHALO Mass Function Models 



mo 

(Me) 



2.0 
1.9 
1.8 
2.0 
2.0 



IQ-^ 2.5 X IQi'^ 
10-6 9.2 X IQi"' 
10"'^ 3.3 X IQi^ 
1 2.5 X IQif 
2.5 X 10^2 



Mtot 
(Me) 



10 



9.3 X 10" 

3.2 X 10" 
2.1 X 10" 
5.8 X 10" 

1.3 X 10^2 



(Me) 



/u 



0.53 
0.19 
0.12 
0.33 
0.73 



7.0 X 10" 0.40 

1.2 X 10" 0.070 

3.3 X 10^" 0.018 
3.5 X 10" 0.20 
1.0 X 10^2 0.60 



Note. — The total number (A^tot), mass (Mtot), and mass 
fraction (/tot = Mtot/M^^^i^) by extrapolation of the subhalo mass 
function with slope a and cutoff mo, normalized to give / = 0.1 
in the interval 10"^ < M/Afjjaio < 10"^. Mu and fu are the mass 
and mass fraction of all subhalos below VL-II's resolution limit of 
~ 10^ Me. 

explanation of our calculation of this background is in- 
cluded in the Appendix \X[ 

We find that our use of an anti-biased radial distribu- 
tion leads to a diffuse subhalo flux that is almost inde- 
pendent of the viewing direction. The median galacto- 
centric distance of a subhalo (i.e. the radius enclosing 
half of all subhalos) is about 200 kpc in the anti-biased 
ca se, but only 100 kp c for the unbiased distribution used 
bv lPieri et"al] ()2008h . The fraction of subhalos within 8 
kpc (within r^L-" = 21 kpc) is 7 x 10"'' (0.01) in the 
anti-biased case, and 0.02 (0.1) for the unbiased distribu- 
tion. In the unbiased case, subhalos within 8 kpc of the 
Galactic center contribute about 90% of the subhalo dif- 
fuse flux towards the Galactic center, whereas they only 
make up 40% of the flux in the anti-biased case. The 
shift towards larger distances also leads to an overall re- 
duction in the amplitude of the flux. 

The final background considered here is due to annihi- 
lations from the smooth host halo. For this component 
we simply use the angular flux distribution calculated 
from all simulated particles that don't belong to any sub- 
halos. Since higher numerical resolution would have re- 
solved some of this DM mass into individual subhalos, 
whose contribution we have accounted for above, we uni- 
formly reduce the smooth halo flux by a factor (1 — /„). 
where /„ is the mass fraction below 10^ Mq (last column 
in Table El). 

3.2. Particle Physics Parameters 

The particle physics dependence of the annihilation 
signal (Eq. [T|) enters through three factors: M^, the 
mass of the DM particle, (crw), the thermally aver- 
aged velocity-weighted annihilation cross section, and 
dN^/dE, the photon spectrum resulting from a single 
annihilation event. The physical nature of DM is cur- 
rently unknown, and a plethora of particle physics mod- 
els have been proposed to explain its existence. It should 
be noted that not all of these models result in a DM 
particle capable of annihilating, but those models are 
not of interest for the present work. Instead we con- 
sider here only the class of models in which the DM is a 
weakly interacting massive particle (WIMP), such as the 
neutralino in supersymmetric extensions of the standard 
model or Kaluza-Klein excitations of standard model 
fields in models with universal extra di mensions (for a 
recent review of particle DM theories, see iBertone et al.l 
l2005aD . 

For any given class of model it is possible to deter- 
mine a range of Af^ and (cru) resulting in a current relic 
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Fig. 6. — The gamma-ray spectrum dN/dE resulting from 
WIMP annihilations into bb quark pairs, calculated with DarkSUSY 
IIGondolo et al.ll2004l) . for between 50 and 500 GeV iri incre - 
ments of 50 GeV. The dashed line shows the lFornengo et al.l (120041 ) 
fitting function for = 100 GeV and the dotted line the expected 
effective area of the GLAST/LAT detector (right abscissa). 



DM density that is consistent wit h the WMAP mea- 



surement of ^l^h^ 



0.1105^°;^°^^ (|Spergel et al .l [2b07 1. 



Typical values for are 50 GeV up to ~ 1 TeV, and 
a simple estimate of the cross section is (av) = 3 x 
lO'^ cm^s'VSly/i^ w 3 X 10~2'5 cm^ s~^ ()Jungman et all 
'199(3'). However, this naive relation can fail badly 
(Profumo 2005,), and a much wider range of cross sec- 



-24 



for < 200 GeV 



tions, up to ( av) ~ 10~"" cm" s 

(e.g. Fig. 17 in lColafrancesco et al.|[2006f ). should be con 
sidered viable. In this work we consider values of 
from 50 to 500 GeV, and (av) from 10"^^ to 10"^^ cm^ 

WIMP DM particles can annihilate into a range of 
different particle pairs, including quarks, leptons (e.g. 
r's). gauge bosons {Z^, W^), gluons, and Higgs par- 
ticles. The subsequent decay and hadronization of the 
annihilation products leads to a spectrum of gamma- 
rays, mostly resulting from the decay of tt^ mesons. 
The shape of this spectrum is almost independent of 
the annihilation branch (with the exce ption of the r 
branch, which has a har der spectrum, iCesarini et al.l 
12004 iFornengo et al.ll200l . and can be calculated using 
Monte-C arlo simulation, for e xample with the PYTHIA 
package (jSiostrand et al|[200ll ). Above x = [E^/M^) = 
0.01 the spectrum is typically well fit by a function of 
the form aQX^ ^-^ exp(—aix), with (ao,ai) depend ent on 
the final state ([Bergstrom et a l. 1998: Koushiap pa s et al.l 
[200l iFornengo et al.ll2004f ). GLAST has non-negligible 
sensitivity down to ~ O.lGeV {x < 0.01), where the 
spectrum has turned over and the x^^^ dependence is no 
longer a good match. Here we calculate dN/dE assum- 
ing a 100% branching ratio into bb quarks (jBaltz et al.l 
l2C)07t iPieri et all I2008D and calcu late the spectrum d i- 
rectly using the DarkSUSY code (|Gondolo et al.ll2004D . 
Annihilation into two photons or a photon and a Z'^ , re- 
sulting in a monochromatic gamma-r ay line signal, is also 
possi ble, but is one-loop suppressed |Bergstrom fc Ullicl 
Il997f ) and we do not include such channels here. Fig. [6] 
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shows the bb spec trum for My between 5 and 200 GeV, 
together with the iFornengo et alJ (|2004[ ) fitting function 
and the GLAST/LAT A^s{E). For = 100 GeV, this 
spectrum results in 32.6, 13.4, and 4.58 gamma-rays per 
annihilatio n above 0.1, 1.0, a n d 3.0 GeV, respectively. 

Recently iBringmann et all (|2008f l presented new cal- 
culations of the gamma-ray spectrum from DM annihi- 
lation, including electromagnetic radiative corrections to 
all leading annihilation processes. In some cases the in- 
clusion of internal bremsstrahlung (IB) photons leads to 
large enhancements of the gamma-ray fluxes, and also 
sharpens the spectral feature at the mass of the DM par- 
ticle. Above energies of 0.6M^, IB photons typically in- 
crease the gamma-ray flux by no more than a factor of 
two over the range of considered here. For the much 
lower energy threshold of 3 GeV employed in this study, 
the enhancement from IB photons is negligible, and thus 
we have not included IB photons here. 

3.3. EGRET Constraints 

Given the theoretically motivated modifications and 
additions to the simulated flux distribution discussed in 
the previous sections, we must ask whether our models 
still satisfy existing constraints from EGRET on the dif- 
fuse gamma-ray flux as well as limits on a Galactic center 
gamma-ray point source. 

Integrating Eq. [13] from 30 MeV to 10 GeV, we obtain 
a specific intensity of 5.91 x 10~^ cm^^ s^^ sr^^ for the 
extra-galactic background in the EGRET energy range. 
Over this same energy range, DM annihilations result in 
30-50 gamma rays per event for = 50 — 500 GcV(see 
Section . The diffuse flux plotted in Fig. [5] can thus 
be converted to units of cm~^ s~^ sr^^ by multiplying 
by (3.09 X lO^icm/kpc) x (30 - 50) x {av)/M^. Even 
the strongest diffuse background model considered here, 
(a,mo/M0) = (2.0,10-12)^ remains below the EGRET 
extragalactic background for (cru)/A/^ < 10~^* cm^ s~^ 
GeV~^, i.e. over the full range of and {<tv) that we 
c onsider in the following a nalysis . 

iHooper fc Dingu^ (j2004[ ) performed a re-analysis of the 
EGRET Galactic center data with an improved energy 
dependent point spread function and found no evidence 
of a point source at the location of Sag A*. For 
between 50 and 500 GeV they determined a 95% con- 
fidence upper limit of '-^ 10~* cm~^ for gamma 
rays above 1 GeV. The total flux from the innermost 
30 arcmin (EGRET 's angular resolution) from a cen- 
tral cusp matched to VL-IFs density profile at rconv 
is equal to 5.7 x 10^^ GeV^ kpc cm~^ for a slope of 

— 1.0. In the energy range from 1 to 30 GeV, DM an- 
nihilations produce 8 — 33 gamma rays per event for 

= 50 — 500 GeV, and the central flux exceeds the 
EGRET limit if (av) > 4.1 x 10-26(Mx/100 GeY)^ '^^ 
cm^ s~^. Most of the parameter space we considered here 
satisfies this limit, but models with < 80 GeV and 
(av) = 3 X 10~^^ cm'^ s~^ violate it. For an inner slope of 

— 1.2 the constrait is stronger, and the EGRET limit is 
already exceeded if (av) > 1.5 x 10-26(Mx/100 GeY)^-^^ 
cm^ s~^. In this case a DM particle mass greater than 
170 GeV (for (av) = 3 x 10"^^ cm^ s'^) would be re- 
quired in order to remain below the EGRET point source 
limit. For the uncorrected central flux, which may in fact 
be closer to reality if dynamical processes have removed 



a significant amount of DM from the Galactic center, the 
limit is relaxed to (av) > 1.5 x 10-'^'^{M^/100 GcY)'^-^'^ 
cm^ s~^, in which case all models considered in this work 
would satisfy the limit. 

3.4. Detectable subhalos 

We have now assembled all necessary ingredients 
to convert our simulated particle distribution into 
a quantitative prediction for the number of subha- 
los detectable with GLAST. To recap: Ng, the ex- 
pected number of gamma-rays from a subhalo in 
one GLAST 9 arcmin "pixel" , is calculated by sum- 
ming B{M)pimp/ (Andf) over all the subhalo's par- 
ticles falling within the pixel and multiplying by 
Tcxp(crw)/(2M2) J^'^ dN-y/dEAcs{E)dE. The total sig- 
nal from a subhalo is obtained by summing Ns over all 
S/N > 1 pixels covered by the subhalo. This flux is then 
compared to ^/^N^,, the square root of the expected 
number of background photons, made up of an isotropic 
extragalactic astrophysical component (as measured by 
EGRET), an astrophysical Galactic background (calcu- 
lated with GALPRO'P), and two diffuse DM annihila- 
tion backgrounds, one from undetectable subhalos and 
one from the smooth host halo. The resulting detection 
significance S = Y^ ^s/y/Yl depends critically on the 
mass of the DM particle (smaller Af^ is better) and on its 
annihilation cross section (larger (av) is better) . At fixed 
and (av) , S depends on the normalization, the slope 
a, and the low mass cut-off mo of the subhalo mass func- 
tion through the calculation of the boost factor and the 
two dilfuse annihilation backgrounds. Note that changes 
in a and toq affect S in opposite ways: increasing the 
abundance of subhalos below VL-II's resolution limit, 
by increasing a or lowering toqi raises the boost factor 
(good for detection) as well as the background from un- 
detectable subhalos (bad for detection). 

Due to departures from spherical symme try in the host 
halo mass and subhalo spatial distribution (|Kuhlen et al.l 
|2007| ). the number of detectable subhalos will depend 
on the observer position. For an observer located at 
8 kpc along the major axis of the host halo ellipsoid 
the diffuse host halo background will be higher than 
for an o bserver at the same distance a long the mi- 
nor axis (jCalcaneo-Roldan fc Mo"ot3 |2000D . Individual 
bright subhalos may also become undetectable from cer- 
tain vantage points, if they fall behind the Galactic disk 
or center for example. For these reasons we have per- 
formed our analysis for ten randomly drawn observer po- 
sitions, in addition to our fiducial observer located along 
the host halo intermediate axis. Our results are pre- 
sented in Fig. [71 where we plot N^, the number of VL- 
II subhalos exceeding 5 = 5, as a function of and 
(o'u)_26 = (o'u)/(10~^^ cm'^ s^^), for different choices of 
a and ttiq. 

We find that ranges from a few up to multiple 
tens over the range of and (av) we have consid- 
ered. In our fiducial model, (M-^/GeV, (cti')_26, a, mo) ~ 
(100, 3, 2.0, 10^^), an observer positioned along the host 
halo's intermediate axis would be able to detect N5 = 19 
subhalos. For the set of ten randomly placed observers 
we find a range of iVs from 13 to 19. Lowering mo 
fro m IQ-^Mp^ to IO^^Mq, the smallest value found 
by iProfumo et al.l ()2006[ ). results in a factor ^ 2 in- 
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Fig. 7. — N^s, the number of simulated subhalos exceeding cS = 5, as a function of the DM particle mass for (av) = 3 X 10~^^ cm^ 
s^^ (top) and the cross section {av) for = 100 GeV (bottom). Dependence on the subhalo mass function cutoff mass mo for slope 
a = 2.0 (left) and on a for mg = 10~® Mq [right). The a = 1.8 case is almost identical to o = 1.9 and has been omitted from this plot. 
The shaded regions indicate the range of TVs for ten randomly chosen observer locations and the solid lines refer to an observer placed along 
the intermediate axis of the host halo ellipsoid. The dotted lino is the case without a boost factor. 



crease in N^. In this case, assuming a = 2.0, more 
than ten subhalos should be detectable even with 
as high as 300 GeV (for {av)-26 = 3) or (av) as low 
as 10-26 cm3 s-i (for = 100 GcV). If, on the 
other hand, no subhalos exist with masses less than 
IMq, then iVg drops below 10 at sa 160 GeV (for 
{av)-26 = 3) or {av)-26 ~ 2 (for = 100 GcV). 
Reducing a below the critical value of 2.0 suppresses 
A^5 by about a factor of 2/3. At a fiducial value of 
(M^, {av)-26) = (100, 3) we find iVg = 13 - 19, 26 - 40, 
11 - 14, 10 - 14, and 9-14 for (a, mo) = (2.0,10-6), 
(2.0,10-12)^ (2.0,10°), (1.9,10-6), and (1.8,10-6), re- 
spectively. In the most favorable case considered here, 
(M^/GeV, (c^^>)_26,a,"^o) = (100,10,2.0,10-^2), the 
number of detectable subhalos exceeds 100. The de- 



pendence of A^5 on niQ and a is primarily driven by 
their effect on the boost factor: lowering the abundance 
of DM subhalos below VL-II's resolution limit reduces 
the boost factor and lowers N5. For comparison we 
have also plotted results without a boost factor for the 
{a, mo) = (2.0, 10-^) case. This case is similar to the 
mo — 1 Mq case, and even without invoking such a sub- 
substructure boost factor multiple subhalos should be 
bright enough for detection for most of the parameter 
space probed here. 

The set of models we have considered inc ludes a model 
simila r to the one recently proposed by iHooper et al.l 
()2007f ) to explain the excess microwave emission around 
the Galactic c enter measured b y WMAP (so - called 
"WMAP haze" . iFhlkbeinal [200I . iHooper et all (pOOTi ) 
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Fig. 8. — The number of detectable (<S = 5) subhalos with more 
than A^pix detectable pixels versus Afpjx, for three different choices 
of M-^ (assuming (av) = 3 X 10~^^ cm^ The shaded regions 

show the range of A'^(> A'^pix) for ten randomly chosen observer 
locations and the solid lines refer to an observer placed along the 
intermediate axis of the host halo ellipsoid. 
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Tp [degrees] 

Fig. 9. — The cumulative angular distribution of detectable (5 > 
5) subhalos, for three different choices of (assuming (av) = 3 X 
10~^^ cm'^ i/i is the angle from the direction of the host center 

to the brightest pixel in each subhalo. The shaded regions show the 
range of N{> tp) for ten randomly chosen observer locations and 
the solid lines refer to an observer placed along the intermediate 
axis of the host halo ellipsoid. The thin solid line represents an 
isotropic distribution normalized to the fiducial observer's location. 



showed that the intensity and angular distribution of 
the WMAP haze can be modeled as synchrotron emis- 
sion from relativistic electrons and positrons generated 
in dark matter annihilations, assuming a cusped cen- 
tral DM density profile of slope —1.2, a DM particle 
mass of 100 GeV, and an annihilation cross section 
of {av) = 3 X 10~^^ cm^ s~^. If this explanation for 
the WMAP haze turns out to be correct, we predict that 
GLAST should be able to detect about 15 subhalos. 

In Fig. [5] we present the cumulative size distributions of 
aU detectable subhalos, for = 50, 100, and 200 GeV 



at fixed {av) — 3x 10~^^ cm'^ s~^. The size here is defined 
as the number of S/N > 1 pixels (A^pix) contributing to 
the total subhalo signal exceeding the detection thresh- 
old of iS = 5. Recall that we have matched the pixel size 
to GLAST's angular resolution of 9 arcmin, such that 
each pixel corresponds to a solid angle of 4.363 x lO^^sr. 
Again we plot distributions for the fiducial observer on 
the host halo's intermediate axis as well as the minimum 
and maximum over all ten randomly drawn observer po- 
sitions. Detectable subhalos are typically extended, with 
a median of TVpix of 13. Only about 5% are detected in 
only one pixel. This should aid in discriminating between 
subhalos lit up by DM annihilation and conventional as- 
trophysical sources (e.g. gam ma-ray pulsars) w hich will 
often appear as point sources (jBaltz et al.ll2007l ). 

The angular distribution of detectable subhalos is 
shown in Fig. [SI where we plot the cumulative num- 
ber of detectable subhalos more than tfj degrees from the 
direction towards the host center. For the fiducial ob- 
server we have overplotted an isotropic distribution (thin 
solid line), in which detectable subhalos are uniformly 
distributed over the whole sky. The actual distribution 
indicates a slight excess at large angles compared to an 
isotropic distribution. However, a KS test shows that 
this discrepancy is not statistically significant. For only 
two of all eleven observer positions can the null hypoth- 
esis of an isotropic distribution be rejected at more than 
90% confidence, and we conclude that the distribution of 
detectable sources on the sky is consistent with isotropy. 

Next we consider the differential distribution of the 
detection significance S, Knax, and the distance to the 
observer, of the subset of detectable subhalos, see Fig.fTOl 
In order to obtain better statistics for this analysis we 
chose a slightly higher cross section of {av) = 5 x 10^^^ 
cm^ s~^, resulting in 51, 29, and 12 detectable subhalos 
for = 50, 100, 200 GeV, respectively. The histograms 
represent the results for the fiducial observer position 
and the error bars indicate the range of values for the 10 
randomly drawn positions. 

The S distribution is slightly peaked towards lower S 
values, with about two thirds of all detectable subhalos 
having S < 10. About 10% of subhalos should be de- 
tectable at very high significances, with S exceeding 25. 
For these, the median Knax is 24 km s~^, indicating that 
the most detectable subhalos tend to be the more mas- 
sive ones. These would be m ore likely to host a dwarf 
galaxy (e.g. lMadau et al.ll2008f ). which m some cases may 
have so far eluded discovery due to their ultra-faint op- 
tical surface brightness. The Vmax distribution, however, 
reveals that a siginificant number of lower VJnax subhalos 
should also be detectable. From 5 to 25 kms~^ the T4iax 
distribution is approximately flat, and in the = 50 
GeV case it even exhibits a pronounced peak towards 
lower Knax- No detectable subhalos are found with Vn^ax 
less than 5 km s~^, which is probably due to a lack of 
numerical resolution in these low mass halos. It appears 
likely that the total number of detectable subhalos has 
not yet converged in our simulations. We explore this 
further in Section [531 below. 

The distance distribution of detectable subhalos is 
broadly peaked around the median of 40 kpc from the 
observer. For the M^^ = 100 GeV case, it is consis- 
tent with a simple log-normal distribution, centered at 
32 kpc with a width of aiog-^„D = 0.45. About 80% of all 
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Fig. 10. — Differential distribution of deteetion significanee S 
(top), Vmax (center), and distance to tlic observer (bottom) of tfie 
subset of detectable (S > 5) subhalos, for three different cfioices of 
Mx (assuming (crv) = 5 X 10~^^ cm^ s~^). The solid lines are for 
the fiducial observer placed along the intermediate axis of the host 
halo ellipsoid and the error bars indicate the range of values for 10 
randomly chosen observer locations. For the distance distribution 
the best-fit Gaussian has been overplotted for the M-^ = 100 GeV 
case. 



detectable subhalos have distances between 10 and 100 
kpc, and only 7.5% are closer than 10 kpc. 

3.5. Convergence? 

In Section [Ol wc discussed the diffuse background re- 
sulting from undetectable subhalos, but we should also 
consider whether any subhalos below VL-II's resolution 
limit of 10^ Mq might be detectable as individual ex- 
tended sources. The fact that the Vmax distribution of 
detectable subhalos remains flat (or even rises slightly) 
down to the smallest visible subhalos hints that the num- 
ber of detectable subhalos has not yet converged in our 
simulations. Further evidence comes from directly com- 
paring the VL-II results with the lower resolution VL-I 
simulation. In VL-I we find far fewer detectable subha- 
los. In fact, with (av) = 3 x 10~^^ cm"^ s~^, at most 
15 subhalos are visible, and that only for the most opti- 
mistic assumptions. Typically only VL-I's most massive 



subhalo reaches 5 = 5. The likely explanation for the 
dearth of detectable subhalos in VL-I is that the lower 
numerical resolution resulted in a significant suppression 
of substructure close to the center. In VL-I only 42 sub- 
halos with at least 200 particles are closer than 50 kpc 
from the halo center, whereas in VL-II this number is 
362, of which 108 are more massive than 4.2 x 10^ Mq 
(200 VL-I particles). Subhalos in VL-II also typically 
have higher central densities: the median of the maxi- 
mum density in all subhalos with Mgub > 4.2 x 10^ Mq 
is 5.4 times higher in VL-II than in VL-I. This will lead 
to higher fluxes and hence increased detectability. 

These considerations motivate a calculation of the to- 
tal expected number of detectable subhalos below the 
simulation's resolution limit, including the full hierar- 
chy of substructure all the way down to the free stream- 
ing cutoff mo- Microhalos, the remnants of the earli- 
est collapsed DM structures, are of particular interest 
here, as they have recently been suggested as possible 
GLAST sou rces that may even exhibit proper motion 
(jKoushia p pas 2006 1 . In the following discussion we focus 
on resolved subhalos, because above the VL-II resolution 
the unresolved component is small, and, as we now show, 
the fraction of unresolved, detectable sources will further 
decrease with mass. 

In the previous section we showed that in our sim- 
ulations 95% of all detectable subhalos are extended 
sources covering more than one pixel. The fraction of 
point-like source is likely to become even smaller for less 
massive systems. For an NFW density profile. 87.5% 
of the annihilation luminosity is generated within the 
scale radius, hence rg/D is a good proxy for a subhalo's 
angular size. Assuming an angular resolution of 9 ar- 
cmin, the maximum distance out to which a subhalo 
will a ppear extended is Umax = ''s/9 arcmin = 380 r,,. 
For a iBullock et"an (|200l concentration-mass relation 
c(M) cx M-°-°'^^ for halos with masses below 10'' Mg. 
In this case 

£»max oc r, cx R/c CX M^^^ M^'"^^ cx A/°-358. (15) 

Specifically, assuming that local subhalos are three times 
more concentrated than field halos (Eq. [Ti|) , we find 



Dr. 



0.6 



M 



10-6 M,. 



0.358 



pc. 



(16) 



At distances exceeding I?max, sources appear point-like 
and the surface brightness in one pixel simply scales as 
L/D^. For microhalos with M < IMq, L cx M°-^^ 
and hence the maximum surface brightness of such un- 



resolved sources scales with mass as L/D. 



2 

max 



93 



(•j^^o.358j-2 ^ ^0.214^ resulting in a smaller fraction of 
detectable point-like sources for lower mass subhalos. 

The differential mass function of extended (i.e. 
GLAST-resolvable) subhalos scales approximately as 
diVext/dM cx dnsn^,/dMDl^^ ^ implying 
that even for a = 2.0 the total number TVcxt will be domi- 
nated by the most massive subhalos, i.e. those just below 
VL-II's resolution limit. The results of a more realistic 
calculation (details in Appendix IB|) . taking into account 
the radial concentration bias and a log-normal concen- 
tration distribution at a given mass, are consistent with 
this simple scaling and arc shown with the thick upper 
lines in Fig. [11] Integrating this differential extended 
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Fig. 11. — The differential mass function of extended (i.e. 
GLAST- resolvable, Ts/D > 9 arcmin, thick upper lines) and ob- 
servable (ts/D > 9 arcmin and "I" > 1 GeV^ kpc cm~^ sr~^, thin 
lower lines) subhalos below VL-II's mass resolution of 10^ Mq, for 
« = 2.0 (solid), 1.9 (dotted), and 1.8 (dashed). 

source mass function from mo = 10^^ Mq to 10^ Mq, 
we find iVe^t = 1090,324,115 for a =2.0, 1.9, and 1.8, 
respectively. 

Of course not all of these extended subhalos would ac- 
tually be detectable by GLAST. We can extend the above 
calculation by imposing an additional constraint on the 
subhalo surface brightness. For an NFW density profile 
the surface brightness within Ts is given by 

L _ L _ 0.875 2 

Since p., cx cV/(c), ^ oc Rc^/f{cy cx M^/^e'/f{cf. 
Thus for a halo of a fixed mass, the surface brightness 
has a very steep concentration dependence, $ cx c'^-3i6 
for M < 10^ Mq. For this reason it is very important 
to account for the radial bias and the log-normal scatter 
in subhalo concentrations. A conservative choice for the 
limiting surface brightness for detection is $o = 1 GeV^ 
kpc cm^^ sr~^, which is a bit below the minimum cen- 
tral brightness of all detectable subhalos in our study. 
The thin lower lines in Fig. [TT] show how this additional 
constraint lowers the expected number of detectable ex- 
tended microhalos. Lower mass halos are more strongly 
affected by a surface brightness cut, since $ scales as 
^^0.225 fQj. ^ M-°-°25. In total we find the surface 
brightness constraint lowers the expected number of de- 
tectable subhalos with masses below 10^ Mq by about 
an order of magnitude to 82, 32, and 13 subhalos for 
a =2.0, 1.9, and 1.8, respectively. Thus we conclude 
that subhalos below VL-IFs resolution limit contribute 
at most 3/4 of the expected total number of detectable 

extended sources. 

Contrary to previous results (jKoushiappasI |2006[ ) . we 
find that GLAST is unlikely to be able to detect and 
resolve individual microhalos, defined here as subhalos 
with masses less than IMq. According to our cal- 
culations, only 0.40, 0.042, and 0.0044 subhalos with 



masses between 10^^ and 1 Mq would satisfy both an- 
gular size and surf ace brightness const raints. The dis- 
agreement between iKoushiappasI (|2006[ ) and our results 
is explained by our lower local subhalo abundance = 
{A'Pdn/dM)/phostiRQ) = 8 X 10"'' versus his 0.002), 
higher local concentrations resulting in smaller angular 
sizes, and a more realistic treatment of the expected 
backgrounds, including a contribution from undetectable 
subhalos. For comparison, we find that the local abun- 
dance of microhalos is (10.7, 0.388, 0.0137) pc'^ for 
a = (2.0, 1.9, 1.8) for the anti-biased radial subhalo dis- 
tribution and (189, 6.83, 0.241) pc~'^ in the unbiased 
case. 

4. DISCUSSION AND CONCLUSION 

In this work we have considered the possibility of di- 
rectly observing Galactic DM substructure through the 
detection of gamma rays originating in the annihilation 
of DM particles in the centers of subhalos. Based on a 
hybrid approach, making use of both the highest reso- 
lution numerical simulations available as well as analyti- 
cal models for the extrapolation beyond this simulation's 
resolution limit, we have shown that for reasonable val- 
ues of the DM particle physics and subhalo mass function 
parameters, future gamma-ray observatories, such as the 
soon to be launched GLAST satellite, may very well be 
able to detect a few, even up to several dozen, such sub- 
halos. 

An important overall systematic uncertainty in our 
study is the nature of the DM particle. Our results are 
valid for the case of a weakly interacting relic particle, 
such as the lightest supersymmetric particle in supersym- 
metric extensions of the standard model. If the DM par- 
ticle instead turns out to be an axion, for example, then 
DM annihilation would not occur and our results would 
be irrelevant. Even in the WIMP case, the allowed pa- 
rameter space for the mass and annihilation cross section 
of the DM particle spans many orders of magnitude. In- 
dividual subhalos will only be detectable with GLAST, 
if and (av) fall in an "observable window" , roughly 
given by ~ 50 - 500 GeV for (av) - 1 - 10 x lO^^e 

s — 1 
cm'^ s . 

Another set of uncertainties is associated with our use 
of a numerical simulation. First of all, we have only sim- 
ulated one host halo at very hig h resolution . Stat istical 
studies at lower resolution (e.g. I Reed et"al] l2005f ) have 
found about a factor of two halo-to-halo scatter in the 
substructure abundance. Secondly, an uncertainty arises 
from the importance of DM substructure below the res- 
olution limit of our simulation, due to its boosting effect 
on the brightness of individual subhalos, its contribution 
to a diffuse background from undetectable subhalos, and 
as individually detectable sources. Lastly, the lack of a 
baryonic component in our simulation is another source 
of uncertainty. In addition to the effects baryonic physics 
might have on the DM distribution at the Galactic cen- 
ter (see Section [2. 3p . one may worry that the abundance 
or detectability of nearby subhalos could be reduced by 
passages through the Milky Way's stellar disk. We can 
use the subhalo orbits from our simulations to constrain 
how large an effect such disk crossings would have. 

For this e stimate we have used orbits from the VL- I 
simulation (jDiemand et al.ll2007d iKuhlen et al.ll2007l) . 
as we are still in the process of extracting them from the 
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VL-II dataset. We model the Milky Way's stellar disk 
as an exponential disk with scale radius Rd = 3.5 kpc 
and scale height Zd = 350 pc. Only 8.2% of all VL-I 
subhalos ever entered a central disk region delimited by 
four scale lengths, i.e. R < ARd and \z\ < Azd-^ This 
fraction grows to 22% when only subhalos within 50 kpc 
of the halo center at z = are considered. Only about 
one fifth of all disk crossing subhalos do so more than 
once. The median of T4, the velocity component per- 
pendicular to the disk, is quite large, about 400 km s""'^. 
Since the typical disk crossing times are thus a few Myr, 
much less than the internal dynamical time of subhalos 
even within (~ 100 Myr), we can apply the impulse 
approximation to determine how m uch subhalos will be 
heate d by disk passages (see §7.2 of lBinnev fc Tremaind 
Il987( ). The change in the z-component of the velocity 
of a DM particle belonging to a disk crossing subhalo is 
given by 

A^. = ^|.g.WI, (18) 

where \z\ is the height of the particle above the sub- 
halo's midplane and gz(R) is the gravitational field of 
the disk at radius R, and is approximately equal to 
5^(i?) = 27rGI]ocxp[(i?o - R)/Rd]- Here Eq ~ 75 Mq 
pc^^ is the surface density of the disk in the solar neigh- 



borhood R — Rq. Setting z 
A\4(|z|) = 5 kms"^ 



exp 



rvmax, we find 
Rq — R 



Rd 



(19) 

V500pc J V400kms-i i ■ y J 
We can understand how destructive such kicks will be 
to the subhalo by comparing AVz to the local circular 
velocity at height z. We find that for all disk crossings in 
VL-I the median of (AVz/Vc) is 0.4 at z = Vg, and only 
25% have AVz > at r^. These numbers will be even 
smaller closer to the subhalo's center, since AVz oc z 
and Vc oc ^/r. We conclude that only about 5% of all 
subhalos within 50 kpc today might have experienced 
a significant reduction in central density, and hence in 
annihilation luminosity, from disk crossings. 

Given the substantial uncertainties discussed above, 
one may wonder whether it is even sensible to consider 
DM annihilations from subhalos at all. We believe that 
the mere possibility of detecting DM annihilations with 
GLAST offers such an exciting opportunity to directly 
confirm the existence of a DM particle and to learn some- 
thing about its properties, that it warrants theoretical 
investigations such as the present one. 
We conclude by summarizing our main findings: 

• Numerical simulations indicate that DM is far from 
smoothly distributed throughout the Galactic halo. 
Extremely high resolution simulations, such as the 
Via Lactea series, have shown that this dumpiness 
extends into individual subhalos, resulting in sub- 
substructure. Since the annihilation rate is propor- 
tional to the DM density squared, such dumpiness 
leads to a boost by (1 -I- B{M)) of the annihila- 
tion luminosity compared to a smooth density dis- 
tribution. Any determination of the gamma-ray 

^ For simplicity we aligned the disk with the z-axis, but the 
results are insensitive to this choice. 



brightness of individual subhalos must account for 
this boost factor. In the case of simulated subhalos 
only the portion of the sub-substructure hierarchy 
below the simulation's resolution limit should be 
included in the boost factor calculation. Our ana- 
lytical model of the boost factor, which assumes a 
powerlaw subhalo mass function of slope a and low- 
mass cutof f mp , a resolved subh alo mass fraction of 
10%, and a lBuUock et al.l (j2001f ) like concentration- 
mass relation, results in moderate boost factors of 
a few up to 10 for (a,77io) = (2.0, 10~^). Lower- 
ing mo to 10~^^ raises B{M) by a factor of 2. For 
a = 1.9 B{M) barely exceeds unity. 

• The portion of the substructure hierarchy below 
Via Lactea's resolution also contributes to a dif- 
fuse annihilation background from undetectable 
sources. The magnitude of this background de- 
pends on the spatial distribution of subhalos within 
the host halo in addition to the subhalo mass func- 
tion parameters (a, TOo). We have e xtended the 
model developed bv lPieri et al.l (|2008f l to allow for 
an anti-biased subhalo radial distribution and for 
higher subhalo concentrations closer to the host 
halo center. The resulting background tends to 
dominate the diffuse background from the smooth 
host halo at angles greater than ip = 10° — 60°, 
depending on the choice of (a,mo). At high galac- 
tic latitudes this background can dominate the as- 
trophysical backgrounds (both Galactic and extra- 
galactic) and should be accounted for when deter- 
mining subhalo detectability. 

• By comparing the expected number of gamma-ray 
photons from an individual subhalo to the square 
root of the number of expected background pho- 
tons, we can determine a detection significance 
S for each subhalo. S depends strongly on the 
mass My. and annihilation cross section {av) of 
the DM particle as well as on the subhalo mass 
function parameters {a, mo). We find that for 
My ~ 50 - 500 GeV and {av) ~ 1 - 10 x lO'^^ 
cm'^ s~^, a few subhalos, even up to a few tens, 
exceed 5 = 5, and hence should be detectable with 
GLAST. In the most optimistic case we considered 
{My = 100 GeV, {av) = lO^^s cm^ s^^, a = 2.0, 
and mo = 1O~^^M0) almost 100 subhalos would 
be detectable, whereas increasing My to 500 GeV 
or lowering {av) to 10~^^ cm'^ s~^ sharply reduces 
the expected number of sources, especially for low 
a or high mo. 

• For the part ic ular D M annihilation model that 
iHooper et al.l (|2007f l proposed to explain the 
WMAP haze, an excess microwave emission around 
the Galactic center, namely My ~ 100 GeV, {av) = 
3 X 10"^^ cm'^ s~^, and a central host halo density 
cusp of slope -1.2, we find that GLAST should be 
able to detect around 15 subhalos, for reasonable 
choices of (a,TOo)- 

• Th majority of all detectable subhalos are resolved 
with GLAST's expected angular resolution of 9 ar- 
cmin. The median number of pixels is 13, and only 
5% are detected in only one pixel. 
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• For 9 out of the 11 observer locations we consid- 
ered, the angular distribution of detectable subha- 
los is consistent with an isotropic distribution on 
the sky. 

• Those subhalos that are detected with the highest 
significances (5 > 25) tend to be the more massive 
ones, with a median Knax of 24 km s^^. How- 
ever, the Knax distribution of detectable subhalos 
is flat, or slightly rising, towards lower Knaxi down 
to Kiax = 5 km s"^ 

• Detectable subhalos typically lie between 10 and 
100 kpc from the observer. The distance distri- 
bution is consistent with a log-normal distribution 
centered at 32 kpc and a width of aiog-^^o = 0.45. 

• It appears that the number of detectable subhalos 
has not yet converged in our simulations. Higher 
resolution simulations would likely resolve smaller 
clumps, which are more abundant and more likely 
to be found near the observer. We estimate that 
we may be missing up to three quarters of all de- 
tectable subhalos. 

Our results suggest that searching for Galactic DM 
substructure should be an important part of the GLAST 
data analysis in the upcoming years. 
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Fig. 12. — Diagram of the relevant geometry for the calculation of JF^aio (Eg. lASl l. 



APPENDIX 

A. CALCULATION OF THE DIFFUSE FLUX FROM INDIVIDUALLY UNDETECTABLE SUBHALOS 



We follow the methodology presented in lPieri et al.l ()2008( ). and define the diffuse flux from individually undetectable 
subhalos from a solid angle Aft around the direction of observation ip by 

( dcosed(t)[ dM ( d\\^ [ dcP{c-M,R) J'i,^io{M,c,X,An) ^^^^iM,R), (Al) 

JAfl Jm J\.o.s. Jc dM 

where A is the distance along the line of sight and R{4\ A, 6*, (p) is the Galacto-centric distance, given by 

R{il), A, e, 4>) = ^X^ +R%- 2 Ai?0 (cos 6* cos V' - cos sin esinip). (A2) 
The subhalo mass function and spatial distribtution is given by 



dM 



{M, R) = AM- 



(l + i?/rMW)2 o 



(A3) 



and we have normalized it such that the total mass in subhalos with masses between ttiq = 10 A/^aio and mi = 

lO-^Mhalo is 10% of Mhalo, l-C 



dMM / 47ri?2_^(A/, R) dR ^ O.lA/haio- 
Jo dM 



(A4) 



Note that this anti-biased spatial distribution is motivat ed by num e rical simulations (jKuhlen et al.l l2007l : 
(MadaueFaDHQOS;) and differs fr om the distribution chosen by iPieri et all (|2008f ) , who employed an unbiased dis- 
tribution of the form 



dn. 



sub 



dM 



(A/, R) ^ AM' 



: n{R ~ r^injM)) 

(i?/rMW)(l + iJ/^MW)2 



Mp)^ kpc~ 



(A5) 



where a Heaviside function 7i(i? — rniin(Af)) accounts for the tid al destruction of subh alos below a radius r^ia{M). The 
subhalo concentration-mass relation P(c; M, R) is based on the lBullock et all (|2001h model {F = 0.01 and K = 3.75) 
to which we have added a radial dependence of the form 



cl^\M,R) = cr{M) (—), 



0.286 



(A6) 
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which wc have determined directly from our nmnerical simulations (jPiemand et al.ir2005bl . [20081 ). With this scaling, 
subhalos at Rq are three times as concentrated as field halos. Note that we include a log-normal scatter of width 
fiogj,, c ~ 0.14, such that 



P(c; M, R) -— - -y= cxp 

clnlO V2^aiog,„c 



(logiQC-logipc r)^ 
2 o-f 



(A7) 



The fiux in a solid angle Afl around ip originating from a subhalo of mass M at a position (A, 9, (p) is given by 

^,.,o(M, c, A, ^, 0, m = / dcos & dcp' I A- P'i^^^^-Ay.y'.O'A')) (^3) 

Jah Ji.o.s. 47rA'^ 

where the subhalo's density profile is given by an NFW profile with a constant core below Tc = 10^* kpc: 

p(M, C, r) = I ^r,r.)'{+r,r.r > (A9) 

The NFW scale density ps depends only on c, and the scale radius on M and c. We have checked that our results 
are not strongly dependent on the value of Tc- T changes by less than 1% for between 10~^ kpc and 10~^ kpc. 

FigdU shows a diagram of the geometry relevant for the calculation of .Fhaio- Without loss of generality we can 
choose a coordinate system in which </> = for the position p of a given subhalo, and hence J-"haio is independent of </>. 
The distance r from the center of the subhalo at p = (A, 0, 0) to the point q = (A', Q\ 0') within the cone of integration 
is 

r(A, 61, A', e*', <\J) = v/A2 + A'2 - 2AA' (sin Q sin 6" cos 0' + cos B cos 6"). (AlO) 

We are integrating over a very small opening angle Ail = 4.36 x 10~^ sr, corresponding to 0max = 1.2 x 10"'^, and 
hence we are justified in setting sin0 = and cos0 — 1. i.e. placing the subhalo on the line of sight given by -i/i. In 
this case Jiiaio depends only on M, c. A, and Afi, and is given by 

^haio(M, c. A, An) = ^ I dcos 0' / ^ de, (All) 



COS t 



lo (^2 _^ ^,2 _ cos 6") ( 1 + VC^ + ^'^ - 2^^' cos 6*' 
where = A/r^ and = A'/r^. 

B. CALCULATION OF DETECTABLE AND EXTENDED SUBHALO MASS FUNCTIONS 
The mass function of extended subhalos is given by 

dN.es 



/ d cos 9 d(^ 









JdXX^ J dc P{c; M, R) (M, R) H{D„,,^{M, c) - A), (Bl) 



dM 

where A, R, P(c; M, r), and dusuh/dM are defined as in Eq. lAll (with t/' = 0). 7i is the Heaviside step function, and 

1 / M 



£'max = ''s/9 arcmin = 10 kpc— (B2) 

' ^ c{M,R) \1MqJ ^ ^ 

is the maximum distance out to which a subhalo of mass M and concentration c(M, R) would appear extended. An 
additional constraint on the surface brightness $ can be incorporated by introducing a second Heaviside step function, 
7^($-$o). 



